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Abstract 

We have directly measured the ZZj and Z77 couplings by studying pp — ► 
££j+X, (I = e, fi) events at y/s = 1.8 TeV with the D0 detector at the Fermi- 
lab Tevatron Collider. A fit to the transverse energy spectrum of the photon 
in the signal events, based on the data set corresponding to an integrated 
luminosity of 13.9 pb -1 (13.3 pb _1 ) for the electron (muon) channel, yields 
the following 95% confidence level limits on the anomalous CP-conserving 
ZZj couplings: -1.9 < /if < 1.8 (h% = 0), and -0.5 < hf < 0.5 (/if = 
0), for a form- factor scale A = 500 GeV. Limits for the Z77 couplings and 
CP-violating couplings are also discussed. 
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I. INTRODUCTION 



Direct measurement of the ZZ^ and Z77 trilinear gauge boson couplings is possible by 
studying Z7 production in pp collisions at the Tevatron (y/s = 1.8 TeV). In what follows 
these couplings will be addressed to as ZVj, where V = Z, 7. The most general Lorentz and 
gauge invariant ZVj vertex is described by four coupling parameters, h\ , (i = 1...4) [|TJ. 
Combinations of the CP-conserving ( CP- violating) parameters h\ and h\ (h\ and 
correspond to the electric (magnetic) dipole and magnetic (electric) quadrupole transition 
moments of the ZVj vertex. In the Standard Model (SM), all the ZV'j couplings vanish 
at the tree level. Non-zero (i.e. anomalous) values of the hj couplings result in an increase 
of the Z'-f production cross section and change the kinematic distribution of the final state 
particles D. Partial wave unitarity of the general // — > Z'-f process restricts the ZVj cou- 
plings uniquely to their vanishing SM values at asymptotically high energies 0]. Therefore, 
the coupling parameters have to be modified by form-factors hf = hY /(l + s/A 2 ) n , where s 
is the square of the invariant mass of the Z7 system, A is the form-factor scale, and hfo are 
coupling values at the low energy limit (s ~ 0) [[§. Following Ref. || we assume n = 3 for 
hY 3 and n = 4 for h\ A . Such a choice yields the same asymptotic energy behavior for all the 
couplings. Unlike W'j production where the form-factor effects do not play a crucial role, 
the A-dependent effects cannot be ignored in Zj production at Tevatron energies. This is 
due to the higher power of s in the vertex function, a direct consequence of the additional 
Bose-Einstein symmetry of the ZVj vertices 0. 

We present a measurement of the ZVj couplings using pp — > ££j + X (£ — e, //) events 
observed with the D0 detector during the 1992-1993 run, corresponding to an integrated 
luminosity of 13.9 ± 1.7pb~ (13.3 ± 1.6pb~ ) for the electron (muon) data. Similar mea- 
surements were recently performed by CDF [|J and L3 p. 

II. DETECTOR AND EVENT SELECTION 

The D0 detector, described in detail elsewhere 0, consists of three main systems. The 
calorimeter consists of uranium-liquid argon sampling detectors in the central and two end 
cryostats, and provides near-hermetic coverage in pseudorapidity (rj) for |?7| < 4.4. The 
energy resolution of the calorimeter has been measured in beam tests [fJJ to be 15%/y/E for 
electrons and 50%/\/E for isolated pions, where E is in GeV. The calorimeter is read out in 
towers that subtend 0.1 x 0.1 in 77 x (where is the azimuthal angle) and are segmented 
longitudinally into 4 electromagnetic (EM) and 4 — 5 hadronic layers. In the third EM layer, 
at the EM shower maximum, the towers are more finely subdivided, subtending 0.05 x 0.05 
in rj x <p. Central and forward drift chambers are used to identify charged tracks for \rj\ < 3.2. 
The muon system consists of magnetized iron toroids with one inner and two outer layers 
of drift tubes, providing coverage for \n\ < 3.3. The muon momentum resolution for central 
muons (\r]\ < 1.0) is determined to be <5(l/p)/(l/p) = 0.18(p — 2)/p © 0.008p {p in GeV/c), 
using Z — » /i/i events. 

Z^i candidates are selected by searching for events containing two isolated electrons 
(muons) with high transverse energy Et (transverse momentum px), and an isolated photon. 
The ee7 sample is selected from a trigger requiring two isolated EM clusters, each with Et > 
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20 GeV. An electron cluster is required to be within the fiducial region of the calorimeter 
(\r]\ < 1.1 in the central calorimeter (CC), or 1.5 < \r]\ < 2.5 in the end calorimeters 
(EC)). Offline electron identification requirements are: (i) the ratio of the EM energy to 
the total shower energy must be > 0.9; (ii) the lateral and longitudinal shower shape must 
be consistent with an electron shower ||; (iii) the isolation variable of the cluster (/) must 
be < 0.1, where I is defined as / = [E tot (0A) - E E m(0-2)]/£ E m(0-2), £tot(0.4) is the total 
shower energy inside a cone defined by 1Z = J (Arj) 2 + (A(f>) 2 = 0.4, and _E E m(0.2) is the 
EM energy inside a cone of 1Z = 0.2; (iv) at least one of the two electron clusters must have 
a matching track in the drift chambers; and (v) Et > 25 GeV for both electrons. 

The /x/x7 sample is selected from a trigger requiring an EM cluster with Ex > 7 GeV and 
a muon track with > 5 GeV/c. A muon track is required to have \rj\ < 1.0 and must have: 
(i) hits in the inner drift-tube layer; (ii) a good overall track fit; (iii) bend view impact 
parameter < 22 cm; (iv) a matching track in the central drift chambers; and (v) minimum 
energy deposition of 1 GeV in the calorimeter along the muon path. The muon must be 
isolated from a nearby jet (7^_j e t > 0.5). At least one of the muon tracks is required to 
traverse a minimum length of magnetized iron (/ Bdl > 1.9 Tm); it is also required that 
Pt 1 > 15 GeV/c and p T 2 > 8 GeV/c. 

The requirements for photon identification are common to both electron and muon sam- 
ples. We require a photon transverse energy E T > 10 GeV and the same quality cuts as 
those on the electron, except that there must be no track pointing toward the calorimeter 
cluster. Additionally, we require that the separation between a photon and both leptons 
be AlZfy > 0.7. This cut suppresses the contribution of the radiative Z — > ££7 decays [[|. 
The above selection criteria yield four eej and two /i/ry candidates. Figure [I] shows the 
E T distribution for these events. Three ecj and both /i/i7 candidates have a three body 
invariant mass close to that of the Z and low separation between the photon and one of 
the leptons, consistent with the interpretation of these events as radiative Z — > ££ — > ££7 
decays. The remaining candidate in electron channel has a dielectron mass compatible with 
that of the Z and a photon well separated from the leptons, an event topology typical for 
direct Z7 production in which a photon is radiated from one of the interacting partons [[|. 

III. SIGNAL AND BACKGROUNDS 

The estimated background, summarized in Table I, includes contributions from (i) 
Z + jet(s) production where one of the jets fakes a photon or an electron (the latter case 
corresponds to the ee7 signature if additionally one of the electrons from the Z — > ee de- 
cay is not detected in a tracking chamber); (ii) QCD multijet production with jets being 
misidentified as electrons or photons; (iii) ttj production followed by decay of each r to 

We estimate the QCD background from data using the probability, P(jet — > e/7), for a jet 
to be misidentified as an electron/photon. This probability is determined by measuring the 
fraction of non-leading jets in samples of QCD multijet events that pass our photon/electron 
identification cuts, and takes into account a 0.25 ± 0.25 fraction of direct photon events in 
the multijet sample 0]. We find the misidentification probabilities P(jet — > e/7) to be 
~ 10 -3 in the typical Et ranges for the electrons and photons of between 10 and 50 GeV. 
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We find the background from Z+jet(s) and QCD multijet events in the electron channel by 
applying misidentification probabilities to the jet E? spectrum of the inclusive ee + jet(s) 
and e7 + jet(s) data. The background is 0.43 ± 0.06 events. For the muon channel the 
QCD background is estimated by applying the misidentification probability to the inclusive 
fifi + jet(s) spectrum. The estimation of the QCD background from data in the muon case 
also accounts for cosmic ray background. The combined background from QCD multijet 
and cosmic ray events in the muon channel is found to be 0.02 ± 0.01 events. 



The TT7 background is estimated using the ISAJET Monte Carlo event generator fTO 
followed by a full simulation of the D0 detector, resulting in 0.004 ± 0.002 events for eej 
and 0.03 ± 0.01 events for /1/17 channels. 

Subtracting the estimated backgrounds from the observed number of events, the signal 
is 3.57jlf ;gf ± 0.06 for the e&y channel and 1. 95^29 ± 0.01 for the fx/jj channel, where the 
first and dominant uncertainty is due to Poisson statistics, and the second is due to the 
systematic error of the background estimate. 

The acceptance of the D0 detector for the eej and /z/ry final states was studied using 
the leading order event generator of Baur and Berger . It generates 4- vectors for the Z7 
processes as a function of the coupling parameters. The 4- vectors were then processed using a 
fast detector simulation program which takes into account effects of the electromagnetic and 
missing transverse energy resolutions, muon momentum resolution, variations in position 
of the vertex along the beam-axis, and trigger and offline efficiencies. These efficiencies 
are estimated using Z — > ee data for the electron channel. The muon trigger efficiency is 
estimated from the e/x data selected using non-muon triggers. The offline efficiency for the 
muon channel is calculated based on e/x and Z — > fifi samples. The trigger efficiency for ee7 
is 0.98 ± 0.01 while the efficiency of offline dielectron identification is 0.64 ± 0.02 in the CC 
and 0.56 ± 0.03 in the EC. For the muon channel the trigger efficiency is 0.94to;og, arid the 
offline dimuon identification efficiency is 0.54 ± 0.04. The photon efficiency depends slightly 
on Ej, due to the calorimeter cluster shape algorithm and the isolation cut, and accounts for 
loss of the photon due to a random track overlap (which results in misidentification of the 
photon as an electron) and the photon conversion into an e + e~ pair before the outermost 
tracking chamber. The average photon efficiency is 0.53 ±0.05. The geometrical acceptance 
for the electron (muon) channel is 53% (20%) for the SM case and increases slightly for 
non-zero anomalous couplings. The overall efficiency for the electron (muon) channel for 



SM couplings is 0.17 ± 0.02 (0.06 ± 0.01). The MRSD— ' set of structure functions is 
used in the calculations. The uncertainties due to the choice of structure function (6%, as 
determined by variation of the results for different sets) are included in the systematic error 
of the Monte Carlo calculation. The effect of higher order QCD corrections are accounted 
for by multiplying the rates by a constant factor k = 1.34 ||. 

The observed number of events is compared with the SM expectation using the estimated 
efficiency and acceptance. We expect the signal in the e and [i channels for SM couplings to 
be: Sj™ = 2.7±0.3 (sys)±0.3 (lum) and = 2.2±0.4 (sys)±0.3 (lum) events, where the 
first error is due to the uncertainty in the Monte Carlo modelling, and the second reflects 
the uncertainty in the integrated luminosity calculation. The numbers are summarized in 
Table I. Our observed signal agrees within errors with the SM prediction for both channels. 
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IV. LIMITS ON ANOMALOUS COUPLINGS 



To set limits on the anomalous coupling parameters, we fit the observed E? spectrum of 
the photon (Ej.) with the Monte Carlo predictions plus the estimated background, combining 
the information in the spectrum shape and the event rate. The fit is performed for the eey 



and /i/ry samples, using a binned likelihood method fl^ , including constraints to account 
for our understanding of luminosity and efficiency uncertainties. Because the contribution of 
the anomalous couplings is concentrated in the high E^ region, the differential distribution 
da/dEj. is more sensitive to the anomalous couplings than a total cross section (see insert 
in Fig. [I] and Ref. 0). To optimize the sensitivity of the experiment for the low statistics, 
we assume Poisson statistics for each Ej. bin and use the maximum likelihood method to 
fit the experimental data. To exploit the fact that anomalous coupling contributions lead 
to an excess of events at high transverse energy of the photon, a high-.Ej, bin, in which we 
observe no events is explicitly used in the histogram flI2|. The results were cross-checked 



using an unbinned likelihood fit which yields similar results. 

Figure [l] shows the observed Ej, spectrum with the SM prediction plus the estimated 
background for the e + /i combined sample. The 95% confidence level (CL) limit contour for 
the CP-conserving anomalous coupling parameters h§ and hf is shown in Fig. A form- 
factor scale A = 500 GeV is used for the calculations of the experimental limits and partial 
wave unitarity constraints. We obtain the following 95% CL limits for the CP-conserving 
ZZ'y and Z77 couplings (in the assumption that all couplings except one are at the SM 
values, i.e. zeros): 

-1.9 < h§ < 1.8; -0.5 < hf < 0.5 

-1.9 < hj < 1.9; -0.5 < hj < 0.5 
The correlated limits for pairs of couplings (h^ , h\ ) are less stringent due to the strong 
interference between these couplings: 

-3.3 < hf < 3.3; -0.9 < h z m < 0.9 

-3.5 < hl Q < 3.5; -0.9 < hj < 0.9 
Limits on the CP-violating ZV'y couplings are numerically the same as those for the CP- 
conserving couplings. The limits on the /if , hf , and h] Q couplings are currently the most 
stringent. 

Global limits on the anomalous couplings (i.e., limits independent of the values of other 
couplings) are close to the correlated limits for (h^ , h^ Q ) and {h\ Q) h^) pairs, since other 
possible combinations of couplings interfere with each other only at the level of 10%. This 
is illustrated in Fig. |3], which shows the limits for pairs of couplings of the same CP-parity 
(couplings with different CP-parity do not interfere with each other). 

We also study the form-factor scale dependence of the results. The chosen value of the 
scale A = 500 GeV is close to the sensitivity limit of this experiment for the h\ 40 couplings: 
for larger values of the scale partial wave unitarity is violated for certain values of anomalous 
couplings allowed at 95% CL by this measurement. 
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TABLES 





ee7 




Candidates 


4 


2 


Background: 






QCD 


0.43 ± 0.06 


0.02 ±0.01 




0.004 ± 0.002 


0.03 ±0.01 


Total background 


0.43 ± 0.06 


0.05 ±0.01 


Signal 


3.57itgf ± 0.06 


1.95t?;jjjj ±0.01 


SM predictions 


2.7 ±0.3 ±0.3 


2.3 ±0.4 ±0.3 



TABLE I. Summary of signal and backgrounds. 



FIGURES 



-i 




FIG. 1. Transverse energy spectrum of photons in ee^f and fi/i^ events. The shadowed bars 
correspond to data points, the hatched curve represents the total background, and the solid line 
shows the sum of the SM predictions and the background. The insert shows da/dEj, folded with 
the efficiencies for SM and anomalous (/if = 3.0) couplings. 
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FIG. 2. Limits on the correlated CP-conserving anomalous ZZ~f coupling parameters /if and 
hf . The solid ellipses represent 68% and 95% CL exclusion contours. The dashed curve shows 
limits from partial wave unitarity for A = 500 GeV. 
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FIG. 3. Limits on the weakly correlated CP-conserving pairs of anomalous ZV^f couplings: a) 
(/if ,ft-3 ), b) c ) (^foi^2o)j an d d) (^ioi ^40 )• The solid ellipses represent 68% and 95% 

CL exclusion contours. Dashed curves show limits from partial wave unitarity for A = 500 GeV. 
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